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Abstract 

Infrared spectra of Comet 21 P/Giacobini-Zinner (hereafter 21 P/GZ) were obtained using 
NIRSPEC at Keck II on UT 2005 June 03, approximately one month before perihelion, 
that simultaneously measured H 2 0, C 2 H 6 , and CH 3 OH. For H 2 0, the production rate of 
3.8 x 10 28 molecules s ' 1 was consistent with that measured during other apparitions of 
21 P/GZ retrieved from optical, infrared, and mm-wavelength observations. The water 
analysis also provided values for rotational temperature (T rol = 55 +3 /. 2 K) and the 
abundance ratio of ortho- and para-water (3.00±0.15, implying a spin temperature 
exceeding 50 K). Six Q-branches in the v 7 band of C 2 H 6 provided a production rate (5.27 
± 0.90 x 10 25 s' 1 ) that corresponded to an abundance ratio of 0.139 ± 0.024 % relative to 
H 2 0, confirming the previously reported strong depletion of C 2 H 6 from IR observations 
during the 1998 apparition, and in qualitative agreement with the depletion in C 2 known 
from optical studies. For CH 3 OH, we applied our recently published ab initio model for 
the V 3 band to obtain a rotational temperature (48 +10 /-7 K) consistent with that obtained for 
H 2 0. In addition we applied a newly developed empirical model for the CH 3 OH v 2 band, 
and obtained a production rate consistent with that obtained from the V 3 band. 

Combining results from both v 2 and V 3 bands provided a production rate (47.5 ± 4.4 x 
10 25 s' ! ) that corresponded to an abundance ratio of 1.25 ± 0.12 % relative to FI 2 0 in 
21 P/GZ. Our study provides the first measure of primary volatile production rates for 
any Jupiter family comet over multiple apparitions using high resolution IR spectroscopy. 


Keywords: spectroscopy, infrared, comets, composition 



1. Introduction 


Comet 21 P/Giacobini-Zinner (hereafter 21P/GZ) was discovered in 1900 by Michel 
Giacobini, was re-discovered independently by Ernst Zinner in 1913, and was recovered 
in all apparitions since its discovery with the exception of the unfavorable ones in 1907, 
1920, and 1953. Its Tisserand parameter with respect to Jupiter (Tj = 2.46) classifies 
21P/GZ dynamically as a Jupiter Family comet (JFC), the predominant “feeding zone” 
for which is the scattered Kuiper disk, as predicted by the pre-eminent Nice model 
(Gomes et al. 2005, Morbidelli et al. 2008). This comet has the distinction of being the 
first to be visited by a spacecraft - in September 1985, the International Cometary 
Explorer (originally the International Sun-Earth Explorer-3) flew through its plasma tail, 
passing within 8,000 km of the nucleus. Here we report abundances of simultaneously 
observed C 2 H 6 , CH 3 OH, and H 2 O in 21P/GZ during its 2005 apparition, using the cross- 
dispersed high-resolution IR echelle spectrometer (NIRSPEC) at Keck II and applying 
our latest fluorescence models for H?0, C 2 H 6 , and CH 3 OH, including a new empirical 
model for the V 2 band of CH 3 OH. 

Observational studies indicate substantial chemical diversity among comets (cf., Mumma 
and Charnley 2011, and references therein). The most comprehensive work reported 
optical narrow-band photometric observations of radical (photo-dissociation product) 
abundances spanning some three decades (A’Hearn et al. 1995). This survey established 
“typical” and “depleted” classes depending on their abundances of C2 (and C3) relative to 
CN and OH. It was argued this may reflect differences in carbon-chain chemistry (i.e., 
differing abundances of progenitor ices for Co and/or C3), with inferred carbon-chain 
depletion being considerably more pronounced in JFCs compared with nearly isotropic 
comets (NICs; i.e., dynamically-new or long-period returning, from the Oort cloud) (see 
also Fink and Hicks 1996, Schleicher et al. 2007, Fink 2009, Cochran et al. 2012). 

While highly diagnostic, these findings do not necessarily translate to a similar 
dichotomy in native ice abundances, primarily due to uncertainties in the parentages of 
fragmentation products (daughter or granddaughter species). Compared with optical 
observations, parent species have been studied in relatively few JFCs, owing mostly to 
their typically lower overall gas production relative to NICs. However, in recent years 
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advances in technology have permitted more detailed studies of JFCs in both radio (Biver 
et al. 2004; Crovisier et al. 2009) and infrared regimes (Weaver et al. 1999; Mumma et al. 
2000, 2005, 2011; Dello Russo et al. 2007, 2008, 2009, 2011; DiSanti et al. 2007a; 
Kobayashi et al. 2007, 2010; Villanueva et al. 2006, 201 lb; Paganini et al. 2012). 

Through use of modern high-resolution spectrometers, the IR offers the unique capability 
of measuring trace gases simultaneously with H 2 O, the most abundant volatile in comet 
nuclei and the principal driver of cometary activity for heliocentric distances (Rh) less 
than ~ 3 AU. Symmetric hydrocarbons (e.g., C 2 H 2 , CH 4 , C 2 H 6 , etc.) have no electric 
dipole moment and so are not observable at radio wavelengths since pure rotational 
transitions are forbidden, leaving their infrared vibrational bands as the sole means of 
measuring their abundances in comets through remote sensing. 

Comet 21P/GZ is often regarded as the prototypical carbon-chain depleted comet 
(A'Hearn et al. 1995). In this study we investigate whether its observed C 2 depletion 
carries over to primary (parent) volatiles, in particular to C 2 H 6 , which we quantify 
through application of a comprehensive fluorescence model for the V 7 band near 3.35 pm 
(Villanueva et al. 201 la). We present a newly-developed empirical fluorescence model 
for CH 3 OH, incorporating more than 100 lines in the V 2 band that are interspersed with 
C 2 H 6 V 7 emissions. Application of our empirical g-factors for V 2 together with our 
recently published ab initio quantum mechanical fluorescence model for the V 3 band near 
3.52 pm (Villanueva et al. 2012) enables us to retrieve a highly precise production rate 
for methanol in 21P/GZ. 

We compare our results for 21P/GZ with those from the IR (for C 2 H 6 and CH3OH) and 
radio (for CH3OH) during its most recent previous apparition in 1998. Our study 
represents the first such measure of parent volatile abundances at infrared wavelengths 
over multiple apparitions in any JFC. We also compare our abundances with those 
measured in other JFCs and with abundances of C 2 from optical studies. We end with a 
discussion of prospects for future studies of 21P/GZ, beginning with its next favorable 
apparition in 2018. 
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2. Observations and Data Processing 

We obtained spectra of 21P/GZ on UT 2005 June 03.63 using the near-infrared cross- 
dispersed echelle spectrograph (NIRSPEC; McLean et ai. 1998) on the 10-m Keck II 
telescope atop Mauna Kea, HI. NIRSPEC delivers simultaneous coverage of six echelle 
orders (at L-band) at sufficiently high spectral resolving power (XI AX ~ 25,000) using a 
3-pixel (0.43 arc-sec) wide slit to measure emission intensities of individual lines or 
groups of blended lines for a given molecule. One half-night (second half) was granted 
to obtain baseline measurements of Comet 9P/TempeI 1 approximately one month prior 
to the Deep Impact encounter. Comet 21P/GZ was accessible in the east toward the end 
of the night, permitting us to dedicate the final 40 minutes of available clock time to the 
study reported here. Because of this, and also due its modest brightness (m v = 10 - 12), 
only one instrument setting was possible - this was chosen so as to simultaneously 
encompass C 2 H 6 (in order 23), CH 3 OH (orders 22 and 23), and H 2 O (orders 26 and 27; 
see Table 1). 

The telescope was nodded in our standard ABBA sequence, with each ABBA set 
representing four minutes of on-source integration time and with A- and B-beams each 
containing comet signal (the telescope was nodded by ±6 arc-seconds relative to the slit 
center for the two positions). For our observations, the slit was oriented at position angle 
157° (Fig. 1A), approximately perpendicular to the Sun-comet direction (PA = 260°). A 
total of six ABBA sets (24 minutes on source) were obtained. Individual echelle orders 
were processed as described previously (Villanueva et al. 2011a; Bonev 2005; DiSanti et 
al. 2001; Mumma et al. 2001a). 

Signals from A and B beams were combined, and spectra were extracted over 9 rows 
(1.78 arc-seconds) centered on the row containing peak molecular emission intensity in 
each order (Figs. IB- ID; also see Table 2, note ‘b’). The emission line spectrum was 
isolated by subtracting the cometary dust continuum convolved with the modeled 
atmospheric transmittance function (Figs. 2 through 4). Flux calibration was 
accomplished by obtaining spectra of an IR standard star (HR-5487). Absolute line 
fluxes were referred to the top-of-the atmosphere by ratioing measured line intensities to 
the monochromatic transmittance at each Doppler-shifted line-center frequency (based on 
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geocentric velocity A-dot = -5.42 km s' 1 ). Transmittances were modeled using the Line- 
By-Line Radiative Transfer Model (LBLRTM; Clough et al. 2005) and the latest 
HITRAN database with updated parameters (Rothman et al. 2010; see Villanueva et al. 
2011a for detailed discussion and application). We used the NIRSPEC Slit-viewing 
CAMera (SCAM) to keep 21P/GZ in the slit, thereby ensuring that loss of flux (slit 
spillover) was due primarily to seeing. We quantified these losses when correcting our 
nucleus-centered production rates to global production rates (see §4.1). 

3 . Empirical Model for the V2 Band of CH3OH 

3. 1 The Need for a v? Band Model 

Methanol is a molecule that is routinely observed in comets, both at IR and radio 
wavelengths. It has three fundamental vibrational bands in the 3.3 - 3.6 pm spectral 
region, V 3 (parallel; symmetric CH 3 stretching mode, centered near 2844 cm' 1 ), V 2 , and V 9 
(predominantly perpendicular; asymmetric CH 3 stretching modes centered near 2990 and 
2970 cm' 1 , respectively). (The terms “parallel” and “perpendicular” refer, in the 
symmetric top approximation, to the orientation of the transition moment associated with 
the vibrational mode relative to the top axis; Herzberg 1945, p. 414.) Of these three 
bands, V 3 has traditionally been targeted for measuring CH 3 OH in comets, primarily 
because it is relatively free of blends with emissions from other species and also because 
it is spectrally more compact compared with the asymmetric bands, having a single 
strong, well-delineated Q-branch. The V 2 and V9 bands (e.g., see Herzberg 1945, Fig. 
128, p. 425) exhibit several sub-band Q-branches dispersed in frequency and interspersed 
with P- and R-branch lines (see §3.2). 

In addition to providing an alternate means of quantifying CH 3 OH in comets, the ability 
to interpret V 2 band emissions has a practical application. Prior to the commissioning of 
NIRSPEC in 1999, the only available high-resolution IR spectrometer was CSHELL at 
the NASA-IRTF 3-m telescope (Tokunaga et ai. 1990, Greene et al. 1994). This legacy 
instrument revolutionized the field of gas-phase spectroscopy at 1 - 5 pm, albeit with 
relatively small spectral coverage (within a single echelle order) per setting, for example 
only about 7 cm ' 1 in the “cometary organics” region (~ 3.3 - 3.6 pm). CSHELL is 
adequate for measuring the V 3 Q-branch (having a total width of ~ 1 .5 cm' 1 ), however this 
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required a dedicated, stand-alone setting for CH3OH, and a second setting for measuring 
other parent molecules, for example C 2 H 6 through its V 7 band centered near 2985 cm' 1 . A 
fluorescence model for the methanol V 2 band permits prominent emission features of 
C 2 H 6 and CH 3 OH to be measured simultaneously, even for instruments with limited 
spectral grasp such as CSHELL. Although since superseded by NIRSPEC, and more 
recently by CRIRES at VLT (Kaufl et al. 2004), the IRTF retains a unique niche by virtue 
of its daytime capability and its greater flexibility in scheduling observing time compared 
with VLT or Keck, further emphasizing the need for quantifying CH3OH using V 2 
emissions. Regardless of the instrument, adequate interpretation of cometary spectra in 
the region from ~ 2970 - 3000 cm ' 1 requires disentangling emissions arising from 
different molecules, most notably C2H6 and CH3OH, with lesser contributions from CH4 
and OH prompt emission. A major step involves characterizing CH3OH in this region, 
and this is dominated by V 2 band emission. 

3.2 Aspects of our V 2 Empirical Model 

Toward this end, we developed an empirical model that includes 157 individual lines, 
with line assignments and frequencies adopted from jet-cooled laboratory spectra (Xu et 
al. 1997) and lower state energies taken from Mekhtiev et al. (1999). We incorporated 
approximately 70 percent of these (109 lines) to quantify CH 3 OH in order 23 (see Fig. 3 
and related discussion below). Details regarding our adopted notation, and results of our 
model as applied to 21 P/GZ are given in the Appendix. 

Fluorescence efficiencies (g-factors) for our empirical model are based on higher signal- 
to-noise spectra of the Halley Family Comet 8 P/Tuttle, also observed with NIRSPEC. 
This comet exhibited a high abundance ratio CH 3 OH/C 2 H 6 (Bonev et al. 2008; also see 
Bohnhardt et al. 2008, Kobayashi et al. 2010), making it a good “template” for this study. 
Fig. 3A shows the resulting model compared with the 8 P/Tuttle residuals from order 23. 
After accounting for emissions from other species (specifically C2H6, CH4 and OH 
prompt emission), we judiciously identified 18 separate spectral intervals in which 
identified CH 3 OH V 2 emissions are isolated (see Appendix). Values for our empirical g- 
factors were assigned such that the production rate for each interval, based on a 9-row 
extract centered on the nucleus (Q» c ,i»/; see Eq. 2 below), agreed with that obtained 
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independently from application (to 8P/Tuttle) of our ab initio quantum model for the 
simultaneously observed V3 band in order 22 (Villanueva et al. 2012) 1 . These g-factors 
were then applied to the 21P/GZ residuals (after accounting for C 2 H 6 and OH emissions; 2 
Figs. 3B-3C). Results for each spectral interval are presented in Table Al. 

4. Production Rates and Molecular Abundances 

4.1 Methodology for Measuring Production Rates in 21P/GZ 

Global (or total) production rates (Q !0 t) for parent molecules were calculated using our 
well-documented methodology (Villanueva et al. 2011a, Bonev 2005, Mumma et al. 
2003, DiSanti et al. 2001). In applying our formalism, we measure a growth factor (GF), 
relating Q to t for a primary volatile to that based on a 9-row extract centered on the peak 
emission intensity (the “nucleus-centered” production rate, Q nc ): 

Q.« =Q, re GF - 0) 

The growth factor GF (> I) corrects for loss of flux, due primarily to atmospheric seeing 
as explained in §2. In our analysis of 21P/GZ, the spectral regions (intervals) we 
selected for determining our production rates typically consisted of blends of multiple 
lines. Within each chosen interval, Q nc for that interval (denoted Q nc ,,„,) is proportional 
to the measured line flux (F„„, W m' 2 ), divided by the product of rest frequency (v,-, 
cm" 1 ), fluorescence g-factor (gi,/, photon s" 1 molecule' 1 , evaluated at R|, = I AU), and 
atmospheric transmittance (T,) at the Doppler-shifted frequency for each encompassed 
line, summed over all contributing lines (numbering *n*): 3 


1 The production rate for CH3OH initially published for 8P/Tuttle (Q lol = 49.65 x 10 2:> 
molecules s" 1 ; Bonev et al. 2008) was based on Q nc = 33.1 x 10 25 s" ! (using the V3 Q- 
branch intensity) combined with a slit loss correction factor (GF) of 1 .5. The value 
obtained using the new band model for v 3 as presented in Villanueva et al. (2012) 
revises Q nc downward by ~ 12 percent, to 29.0 x 10 25 s' 1 . 

2 Due to its small geocentric Doppler shift at the time of our observations (-5.42 km s" 1 ), 
CH4 was not observable in 21P/G-Z due to severe extinction by corresponding telluric 
absorption. 

3 For the case of only one line within an interval, the quantity in square brackets in Eq. 2 
is equivalent to (F t , n e/T)/ (V] , ne gime.iAu), where Fjj ne /T is the top-of-atmosphere line flux 
as described in §2. 



(2) 


/=i 

In Eq. 2, Ti (s) is the photo-dissociation lifetime at Ru = 1 AU, f(x) is the fraction of all 
molecules in the coma contained in the nucleus-centered beam (assuming native release 
and sphericaliy-symmetric outflow; see Appendix in Hoban et al. 1991), A is expressed in 
meters, and he = 1.99 x I O ' 23 W s cm. The overall Q nc represents the mean of values 
from all intervals, weighted by their respective stochastic errors. For CH 3 OH in order 23, 
the uncertainty in our equivalent empirical g-facior for each interval is also included in 
this weighting (Table Al). 

We calculate our global production rates using GF = 1 .62 (Table 2), based on the spatial 
emission profile summed from several bright water lines (Fig. 1 A). Although noisier, the 
GF for CH 3 OH was in agreement with this value, and that for C 2 H 6 , although somewhat 
larger, was also consistent (within la uncertainty) with that measured for FI 2 O. Details 
are given in Table 2 (in particular, see note d). 

4.2. Results for Individual molecules 
4.2.1. Ii 2 0 

A total of 31 water lines contained within 24 spectral intervals were measured in orders 
26 and 27 (Fig. 2). These sample a sufficiently large range in rotational energy (> 300 
cm' 1 ; see inset to Fig. 2) to allow determination of the rotational temperature with 
extremely high precision (T ro t = 55 +3 /.2 K). Of the encompassed lines, 18 are ortho 
(parallel nuclear spins of the two H-atoms), and 13 are para (anti-parallel spins). 
Analysis of their relative intensities allowed an accurate measure of an ortho-to-para ratio 
(OPR, 3.00 ± 0.15) that is consistent with equilibrated spins, thereby implying a spin 
temperature (T sp i n ) exceeding 50 K (e.g., see Bonev et al. 2007). Our method for 
decoupling T ro i and T S[) j n provides a highly precise (and accurate) measure of the global 
production rate for H 2 0 (Q(H 2 0), 3.800 ± 0,088 x 10 28 molecules s' 1 ). 


Q 


4 M 


nc ,im 


he T, f(x) 
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4.2,2. OH Prompt Emission 

Our observations encompass four strong lines of OH in order 23 (Fig. 3). These lines 
arise through “prompt emission” in the v = 1 - 0 band, radiating within milliseconds from 
vibrationally and rotationally excited states following their production from photo- 
dissociation of H 2 O molecules in the coma (Bonev et al. 2006, Bonev and Mumma 2006, 
Mumma et al. 2001a). Empirical g-factors for more than 20 OH prompt emission lines 
(denoted OH*) from v = 1 - 0 and 2 - 1 bands were derived from NIRSPEC spectra of 
two comets (C/2000 WM 1 LINEAR and C/2004 Q2 Machholz; Bonev et al. 2006). 

However, such a study is problematic for the lines observed here, in that these are 
interspersed with emissions from C 2 H 6 , CH 3 OH, and perhaps other species. In particular, 
our effort to obtain accurate production rates in 21P/GZ, both for C 2 H 6 and for CH 3 OH in 
order 23 is hampered by uncertainties in the g-factor for the OH PI 2.5 2" line at rest 
frequency 2996.892 cm* 1 . This line is blended with C 2 H 6 r Cb, and also with the strong 
A + '"'(K=l->0) Q-branch of the CH 3 OH V 2 band (plus a few weaker V 2 lines; see Table 
Al). There are three reasons for this uncertainty in g-factor: (1) The spectrally confused 
nature of this region precludes isolating the contribution from OH*, and so requires 
interpolation/extrapolation from line g-factors measured from other “quadruplets” such as 
those near 3045 and 2785 cm* 1 ; (2) the intensities of OH* lines depend on the details of 
water excitation and subsequent uni-molecular dissociation (e.g., on the ratio of solar Ly 
a to UV-continuum ratio; see discussion in §4 of Bonev et al. 2006; see also Bonev and 
Mumma 2006 and references therein). Consequently, it is not known how well the value 
for the PI 2.5 2‘ line in 8 P/Tutt!e applies to 21P/GZ, observed approximately 2.5 years 
earlier; and (3) in 21P/GZ, maximum line intensities for C 2 I-L and CH 3 OH were obtained 
from spectral extracts displaced slightly from those for H 2 O (see Table 2, note ‘b'); this 
introduces additional uncertainty in scaling the modeled OH* contribution. 

Nonetheless we estimated their effective g-factors through comparison with OH lines 
from adjacent multiplets. We also estimated the g-factor for PI 2.5 2* empirically from 
NIRSPEC observations of Comet C/1999 S4 (LINEAR) in July 2000. This comet was 
very highly depleted in virtually all primary volatiles, including C 2 H 6 and CH3OH 
(Mumma et al. 2001b), leaving the OH* lines as the only discernable emissions in excess 
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of the continuum in order 23. These two independent approaches led to g-factors for the 
P12.5 2' line in the ratio 1.19 (the one based on LINEAR S4 being larger), exemplifying 
the uncertainties associated with obtaining accurate g-factors for these lines when 
comparing comets observed at different times. Based on this discussion, we exclude this 
region in calculating our “most probable” production rates and mixing ratios for C 2 H 6 
and CH 3 OH; however, including signal in this region does not alter our conclusions 
regarding the abundances of C 2 H 6 and CH 3 OH (see below, and Tables 2 and 3). 

4.2.3. C 2 H 6 

Order 23 is dominated by emissions from ethane (V 7 band) and methanol (primarily the V 2 
band). One motivation for our study of 21P/GZ is to investigate (and, in the case of 
C 2 H 6 , to confirm) whether the depletion of C 2 observed in this comet (A’Hearn et al. 
1995; Fink 2009) carries over to parent molecules. Depleted C 2 H 6 was reported from 
CSHELL observations of 21P/GZ during the 1998 apparition (Weaver et al. 1999; 
Mumma et al. 2000; see §5.1.2). This might be expected if acetylene (C 2 H 2 ) is the 
dominant source of C 2 , as the most efficient means of producing C 2 H 6 is H-atom addition 
to C 2 H 2 on the surfaces of interstellar icy grains at low temperatures (Mumma and 
Charnley 201 1; Rodgers and Charnley 1998; also see § 6 ). 

We obtained production rates from the intensities of seven Q-branches of the C 2 H 6 V 7 
band: r Q 3 , r Q 2 , r Qi, r Qo, p Qi, P Q 2 , 'Xb, and identify these by numbers 1 - 7, respectively, 
in Fig. 3B. Our production rates were based on a full quantum model for the V7 band of 
C 2 H 6 that includes P- and R-branch intensities and contributions from the V7+V4-V4 
combination band (Villanueva et al., 201 la). Because each Q-branch samples a range of 
rotational energies, and because there is considerable overlap in the energies sampled 
among these Q-branches, the resulting production rate is relatively insensitive to T rot . For 
this reason, we adopt T rot = 50 K for C 2 H 6 . This temperature falls between those we 
measured for H 2 O and for CH 3 OH (based on analysis of the V 3 band; see §4.2.4. 1 and 
Fig. 4). 

Global production rates for the seven Q-branches were calculated using their integrated 
line fluxes and applying the growth factor measured for FLO (= 1 .62); these are indicated 
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in Fig. 3D. However, because the g-factor for the blended PI 2.5 2' prompt emission line 
of OH is uncertain (§4.2.2), r Ch is excluded from our quantitative analysis of C 2 H 6 . The 
mean global production rate is then the mean of the remaining six points, weighted by 
their individual stochastic uncertainties (denoted by the error bars in Fig. 3D). We 
obtained a global production rate (Q^Hg), 5.27 ± 0.90 x 10 2 ~ s s’ 1 ) that translates to an 
abundance ratio of (0.139 ± 0.024) x 10' 2 relative to H 2 O (Tables 2 and 3). This is well 
below the current “normal” value (~ 0.6%; Bockelee-Morvan et al. 2005, DiSanti and 
Mumma 2008, Mumma and Charnley 2011), confirming highly depleted C 2 H 6 in 21P/GZ 
and commensurate with its observed C 2 depletion (A’Hearn et al. 1995, Fink 2009; see 
also Table 3). 

4.2.4. CH 3 OH 

4.2.4. 1. Rotational Temperature and Abundance for CTI 3 OH from Order 22 

We identified 36 spectral intervals in order 22 and, applying our new V 3 band model 
(Villanueva et al. 2012), we measured a rotational temperature for CH 3 OH (T rot - 48 +l0 /-7 
K; Fig. 4) consistent with that measured for H 2 O in 21P/GZ. Our production rate for 
CH3OH from order 22 (4.63 ± 0.55 x 10 26 s' ; ) translates to an abundance ratio of (1.22 ± 
0.1 5) x 10 ~ 2 with respect to H 2 O (Table 3). 

4. 2. 4. 2. CH 3 OH in Order 23 - Results using our Empirical Model for the V? Band 
Subtracting modeled C 2 H 6 (and OH*) emissions from the (continuum-subtracted) 
residuals (Fig. 3B) results in a net emission spectrum dominated by methanol (labeled 
“CH 3 OH residuals” in Fig. 3C). Comparison with our empirical model provides an 
additional measure of CH 3 OH production. In Fig. 3E we show the resulting global 
production rates for CH 3 OH in 21P/GZ based on signal contained in each of the 18 
spectral intervals in order 23. The mean of these (excluding interval 5; see §4.2.2) results 
in Q(CH 30 H) = (4.88 ± 0.71) x 10 26 s ’ 1 and therefore an abundance ratio of (1.28 ±0.19) 
x 10' based on CH 3 OH in order 23 alone. (Including interval 5 results in an abundance 
ratio of (1.06 ± 0.17) x 1 O' 2 .) Combining CH 3 OH from orders 22 and 23 (but excluding 
interval 5) results in Q, ol (CH 3 OH) = (47.5 ± 4.4) x 10 25 s' 1 , or (1 .25 ± 0.12) x 10 ' 2 relative 
to H 2 O. 
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5. Abundance Comparisons 

5.7. Comparisons with Other 21P/GZ Observations 

5.1.1. Water Production 

The global production rate we measure for H 2 O in 21P/GZ (3.8, expressed in I0 28 
molecules s' 1 ) compares favorably with the values reported elsewhere. In their optical 
photometric survey, A* Hearn et al. (1995) give Qoh = 3.2 (corresponding to Qmo ~ 3.5) 
at the same heliocentric distance as our 2005 observations (R h =1.12 AU), and A’Hearn 
et al. estimated 4.2 for the maximum Qh 2 o at perihelion (Rh =1 .03 AU). Fink (2009) lists 
a range (1.6 - 4.2) from optical spectra encompassing 1.05 < Rh < 1.21 AU during the 
1985 apparition. Weaver et al. (1999) measured Qh 2 o ~ 2 - 3 from 1 998 for 1 .05 < Rh < 
1 .1 1 AU based on IR spectra using CSHELL. Crovisier et a!. (2002) report a maximum 
QH 20 “ 5.1 at Rh = 1.05 AU, and Crovisier et al. (2009) report Qh 2 o - 3 at Rh - 1.2 AU, 
based on observations of OH from Nan 9 ay during the 1998 apparition. 

5.1.2. Infrared Observations of C 2 H 6 and CH 3 OH in 21P/GZ 

Highly depleted C 2 H 6 in 21P/GZ was reported in two papers from the 1998 apparition 
based on observations with CSHELL. Mumma et al. (2000) reported an abundance ratio 
(C 2 H 6 /H 2 O, 0.22 ±0.13 x 10' 2 ) that is consistent with ours, albeit with larger uncertainty 
due to the lower sensitivity afforded by CSHELL compared with NIRSPEC. Weaver et 
al. (1999) reported a range of (3a) upper limits, the maximum of their range (0.08 x 10' 2 ) 
being only ~ 60 percent of our value (0.14 x 1 O' 2 ; Table 3). Taken together, these two 
papers suggested possible compositional heterogeneity in 21P/GZ. Our observed 
differences in spatial profiles for H 2 O, C 2 H 6 , and CH 3 OH (Fig. 1) may result from 
compositional differences among source regions on the nucleus. (However, given the 
relatively poor signal-to-noise in our profiles for ethane and methanol, and without 
detailed knowledge concerning the release of material, definitive conclusions cannot be 
drawn.) For methanol, Weaver et al. reported a range for the abundance ratio (0.9 - 1 .4 x 
10' 2 ) based on the entire V 3 band intensity and gband = 1.5 x 10 ' 4 photon s" ! molecule' 1 . 
This encompasses our combined value of (1.25 ± 0.12) x 10 ' 2 (from orders 22 and 23), as 
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well as our values from each order separately, and even our value of ( 1.12 ± 0 . 1 1 ) x 10 
when interval 5 in order 23 is included in our two-order mean CH 3 OH abundance. 

5.1.3. Millimeter Observations of CH 3 OH in 21P/GZ 

Biver et al. (2002) included 21P/GZ in a survey of molecular abundances in 24 comets 
based on millimeter wavelength observations with the Institut de Radio Astronomie 
Millimetrique (IRAM) 30-m and other radio telescopes. For CH 3 OH they reported an 
abundance ratio ( 1.6 ± 0.4 x 10 ' 2 ) that implied a range Qcmoii ~ (3.6 - 10 ) x 10 26 
molecules s' 1 , considering their ±lo extremes in abundance along with the range in Qm 2 o 
from the Nan 9 ay observations (3 - 5.1 x 10 molecules s' ; see §5.1.1). Based on this, 
our measured abundance for CH 3 OH falls near the low end of the IRAM 1 o limit. The 
agreement between these results and our measurements is notable, especially considering 
that the radio measurements of CH 3 OH and OH were conducted with different telescopes 
having greatly different beam sizes whose results represent temporal release averaged 
over widely different intervals, whereas our NIRSPEC observations measured both 
CH 3 OH and H 2 O simultaneously and in the same we II -characterized (small) beam. 

5 . 2 . Abundance Comparisons of C2H6 and CH3OH Among Jupiter Family Comets 

An ensemble of comets (both NICs and JFCs) measured with IR spectroscopy reveals an 
emerging taxonomy based on the diversity of primary volatile abundances (specifically 
C 2 H 6 , CH 3 OH, C 2 H 2 , and HCN) among its members. Seven of 14 comets measured 
through 2008 are conditionally classified “organics-normal,” while two are “enriched” 
and two are “depleted,” particularly regarding their abundances of C 2 H 6 and CH 3 OH, for 
which “normal” values are, respectively, approximately 0.6 and 2 percent relative to H 2 O 
(e.g., see Fig. 8 in Mumma and Charnley 2011). A key question regards the extent to 
which comets fall into these three categories as opposed to forming a more continuous 
distribution in volatile organic compositions, and also whether the extreme examples 
(both enriched and depleted) represent end members or even if the “normal” group 
reflects a representative mean composition. Indeed, this classification scheme is being 
challenged as additional comets are characterized. For example, neither 8 P/Tuttle 
(Bohnhardt et al. 2008, Bonev et al. 2008, Kobayashi et al. 2010) nor C/2007 N3 (Lulin) 
(Gibb et al. 2012) fit neatly into any of the three categories. 
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The overall sample is dominated by NICs (delivered from the Oort cloud), however the 
JFC population (primarily from the scattered Kuiper disk) contributes in significant ways. 
An important finding revealed even in this relatively small sample of comets is that 
members from NIC and JFC populations are present both in enriched and depleted 
categories, as predicted by the Nice model (Gomes et al. 2005, Morbidelli et al. 2008). 

Table 3 exemplifies this through abundance comparisons among JFCs. Comet 
17P/Holmes is clearly enriched (Delio Russo et al. 2008) and 73P/Schwassmann- 
Wachmann 3 is depleted (Villanueva el al. 2006; Delio Russo et al. 2007) in both C 2 H 6 
and CH 3 OH, while lOP/Tempel 2 (Paganini et al. 2012) showed mild depletions in both, 
as did 9P/Tempel 1 (except for the material ejected by the Deep Impact event, which 
revealed “normal” C 2 H 6 ; Mumma et al. 2005, DiSanti et al. 2007b). Comet 2IP/GZ was 
highly depleted in C 2 H 6 yet had only mildly depleted (to nearly normal) CH 3 OH, 
6P/d* Arrest (Delio Russo et al. 2009) had mildly depleted C 2 H 6 and normal CH 3 OH, 
while in 103P/Hartley 2 both molecules had abundances consistent with normal (Mumma 
et al. 201 1; Delio Russo et al. 2011; Kawakita et al. 2012). Clearly there are examples 
that do not fit the proposed relatively simple three-tiered classification scheme, even 
based on this small sample of seven JFCs. 

5.3. Comparisons with C 2 Abundances from Optical Observations of JFCs 

An additional comparison involves the potential relation between primary volatile 
abundances (in particular C 2 H 6 for the present study) and C 2 . In an initial optical survey 
(A’Heam et al. 1995), approximately one-third of 85 comets were classified as “carbon- 
chain depleted,” with this classification being considerably more pronounced in JFCs 
compared with NICs. This dichotomy with dynamical class is reinforced by other 
surveys (Fink and Flicks 1996, Fink 2009, Cochran et al. 2012), including an extensive 
taxonomic compilation of well over 100 comets that identifies up to nine compositional 
groupings when abundances of NH, a potentially important tracer for NM 3 , are included 
(Schleicher et al. 2007). 

Five of seven comets in Table 3 have fairly similar values for the ratio C 2 I-VC 2 . One 
exception is 17P/Holmes, however measurements of both C 2 H 6 and C 2 were obtained 
following its spectacular outburst in October 2007, and so uncertainties associated with a 
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rapidly evolving coma might be expected. This leaves 6 P/d’ Arrest as a notable outlier, 
with mild depletion of C 2 H 6 (by approximately a factor of 2 relative to “normal”) but C 2 
near the midpoint of the range for comets classified as “typical” by A’Hearn et ai. (1995). 

While the majority of JFCs in Table 3 exhibit similar C 2 IVC 2 , the inter-relationship 
between these two species is far from secure. C 2 is produced directly from C 2 H 2 with a 
branching ratio of approximately 20 percent (Huebner et al. 1992), but some comets 
produce more C 2 than can be accounted for by C 2 H 2 alone, for example 8 P/Tuttle (Bonev 
et ai. 2008) and C/2007 N3 (Lulin) (Gibb et al. 2012), both of which were depleted in 
C 2 H 2 . An analogous situation may hold for the JFC 6 P/d* Arrest, with strongly depleted 
C 2 H 2 (< 0.052 %, 3a) suggesting the need for an additional source of C 2 , and similarly 
one or more additional sources for CN by virtue of its larger abundance compared with 
that of HCN (Dello Russo et al. 2009). 

In order to reproduce observed profiles of C 2 in Comet 1 P/Hailey, Combi and Fink 
(1997) investigated CHON grains as a possible additional source, and Klavetter and 
A’Hearn (1994) found that ~ 50 % of CN in the coma of 1 P/Halley originated from jets 
of CHON origin. Lederer et al. (2009) proposed sub-micron organic grains (e.g., CHON) 
to explain OH, CN, and C 2 jets in C/1995 01 (Hale-Bopp), and Frey et al. (2005) 
proposed additional sources of C 2 in comets, including photo-degradation of CH 3 CN, 
HC 3 N, and C 2 N 2 . 

Photo-dissociation of C 2 FH can lead to C 2 , however multiple steps are required (including 
production of C 2 H 2 ) and branching ratios are highly uncertain (Helbert et al. 2005, 
Kobayashi et al. 2010). Indeed, Weiler (2012) compared radial column densities of C 2 
and C 3 observed in three comets with a coma chemistry model, and concluded that C 2 H 6 
is not a significant parent of C 2 . The similar values for C 2 I-VC 2 among five of the seven 
JFCs in Table 3 may be coincidental, or could indicate a common source of (and 
efficiency for producing) C 2 and C 2 H 6 (such as C 2 H 2 ) that dominates contributions from 
alternative sources (e.g., CHON) in these five comets. The problem of relating 
abundances of fragment species to potential parents is not resolved, and parentages may 
vary among comets. Clearly, in general both gas and grain sources need to be 
considered. 
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6. Production of C 2 H 6 and CH 3 OH through Grain Surface Chemistry 

Infrared measurements such as those presented here can provide information on 
conditions in the environment where comets formed. Because gas phase formation of 
C2H6 is highly inefficient, a mechanism was proposed (Tielens 1992) involving H-atom 
addition reactions to C2H2 on the surfaces of interstellar grains, and this was later borne 
out by the discovery of abundant C2H6 in comet C/1996 B2 (Hyakutake) (Mumma et al. 
1996) and through further theoretical development (Rodgers and Charnley 1998). The 
now-routine measurement of C2H6 in comets suggests such processing is ubiquitous in 
pre-cometary environments. 

Similarly, gas-phase formation of CH3OH, through reaction of the methyl ion (CH3 + ) 
with water followed by dissociative recombination of the product species (CIHhOI-k*), 
cannot explain observed abundances of CH3OH in interstellar sources (Geppert et al. 
2006; Garrod et al. 2006; Wirstrom et al. 2011). In addition, chemical models that trace 
the evolution of material from dark molecular clouds through proto-planetary disks have 
not found a viable formation pathway for CH3OH in the gas phase (Walsh et al. 2010; 
Willacy & Woods 2009; Willacy 2007). However, in analogy with C2H6, methanol is 
efficiently synthesized through grain surface chemistry in sufficiently cold environments. 
Hydrogenation of CO condensed onto grains, as simulated by laboratory irradiation 
experiments on pure CO ice and CO/H2O ice mixtures at very low temperatures (~ 10 - 
20 K), show this to be an efficient means of producing CH3OH (Hudson and Moore 
1999, Hiraoka et al. 2002, Watanabe et al. 2004). 

7. Implications for Future Studies of 21P/Giacobini-Zinner 

In this context, our measurements of strongly depleted C 2 H6 but only modestly depleted 
CH3OH in 21P/GZ lead naturally to questions regarding its abundances of C2H2 and CO. 
Weaver et al. (1999) reported an abundance of C2H2/H2O < 0.3 - 0.4 % (3a), which does 
not represent a sensitive upper limit given the (current) “normal” abundance of C2H2 
among comets (~ 0.2 - 0.3 %; Mumma and Charnley 2011, DiSanti and Mumma 2008, 
Bockelee-Morvan et al. 2005). Mumma et al. (2000) reported a high abundance ratio 
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CO/H2O (~ 10 %), however this was based on a marginal detection of H2O (below 2a) 

together with an abundance ratio C2H6/CO = 2.1 ± 0.6 %. Our current result for C2H6 

$ 

(0.139 ± 0.024 %) implies an abundance ratio CO/H 2 0 = 6.6 ± 1.1 %, if scaled from the 
Mumma et al. (2000) ratio for CiW(JCQ. Results from these two studies (Mumma et al., 
Weaver et al.), conducted approximately three weeks apart, are at odds. In the latter 
work H 2 O was detected yet neither C2H6 nor CO were seen, and both papers mentioned 
compositional heterogeneity as a possible explanation for these differences. Assuming a 
(modest) Rif 2 dependence, the upper limit reported for Q(CO) at Rh = 2.4 AU from 2005 
Spitzer observations of 21P/GZ (6.0 x 10 26 molecules s' 1 ; Pittichova et al. 2008) scales to 
~ 2.8 x 10 27 molecules s' 1 at R|, - 1.12 AU (corresponding to the observations reported 
here). Based on our measured Q(H20) (3.8 x 10 28 molecules s' 1 ), this corresponds to 
CO/H2O < 7.3 %, consistent with the revised Mumma et al. value (6.6 %). 

Because of its short orbital period (6.6 years), apparitions of 21P/GZ provide widely 
varying observational opportunities. This is particularly true when comparing the 
apparitions in 2012 and 2018. Near perihelion, the geocentric distance A was > 1.8 AU 
in 2012, but will be only ~ 0.4 AU in 2018. Moreover, the geocentric velocity in 2018 
exceeds 10 km s' 1 when A reaches ~ 0.48 AU (and thus is adequate for studying CO and 
CI-I4). Assuming gas production rates (and rotational temperatures) similar to those 
measured in previous apparitions, line intensities in 2018 are expected to be larger by 
factors of three to four compared with those observed in 2005, and by factors of two to 
three compared with those from 1998 (for which A ~ 0.9 - 1.1 AU). Combined with 
anticipated improvements in instrumentation (e.g., higher spectral resolving power and 
sensitivity), the 2018 apparition will enable deep searches for trace species - e.g., 
detection of C2H2 or, if severely depleted, establishing stringent upper limits for its 
abundance - and will allow the question of potential compositional heterogeneity in the 
nucleus to be addressed in considerable detail. 

8. Conclusions 

We used NIRSPEC at Keck II to obtain infrared spectra of 21P/GZ on UT 2005 June 
03.6 that allowed a simultaneous measure of II2O, C2H6, and CH3OH. For H2O, we 
obtained a production rate (3.8 x 10 28 molecules s' 1 ) consistent with that measured during 
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earlier apparitions of 21F/GZ as retrieved from optical, IR, and mm-wavelength 
observations. Our water analysis also provided consistent values for rotational 
temperature (T rot , 55 +3 /. 2 K) and spin temperature (T spm , > 50 K, based on OPR 3.00 ± 
0.15). 

Our results confirm strongly depleted C 2 H 6 and mildly depleted CH 3 OH in 21P/GZ, as 
reported previously from observations during the 1998 apparition. For C 2 H 6 we applied 
our ab intio quantum mechanical fluorescence model for the v 7 band that includes 
contributions from the Q-, P-, and R-branches and also from a combination band. For 
CH 3 OH, we applied our recently published ab initio model for the V 3 band to obtain a 
rotational temperature (T rob 48 +i0 /. 7 K) consistent with the value we measured for H 2 0. 
The first application of our newly developed empirical model for the CH 3 OH v 2 band 
provided a production rate consistent with that obtained from V3. Combining analyses of 
both bands (v 2 and V3) permitted a ~ 10a detection of CH 3 OH in 21P/GZ. Our study 
represents the first measure of primary volatile abundances from ground-based IR 
spectroscopy for any JFC over multiple apparitions. It also provides a comparator for 
more extensive studies of 21P/GZ in future apparitions, beginning with that in 2018. 
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Table 1. Observing Log for Comet 21P/Giacobini-Zinner on UT 2005 June 03 a 


NIRSPEC 

Setting 

KL1 


UT 

Start-End 


Molecule (order) b 


Tint c 

[min] 


14:45-15:24 H 2 0 (26,27) 24 

C 2 H 6 (23) 

CH 3 0H (22, 23) 


a Rh = 1.116 AU, A = 1.454 AU, dA/dt = -5.4 km s* 1 , dRh/dt = -8.8 km S' 1 . Solar phase 
angle ((3) = 44°. 

b Quantified primary molecules and NIRSPEC echelle order in which their emissions 
occur. 

c Total on-source integration time. 
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Table 2. Primary Volatiles in 21P/G-Z on 2005 June 03.6 

Molec. 

Order Band(s) ^ ^ GF^ ^ 

[K] (10 2 5 S'*) (lO^S- 1 ) 

H 2 0 

26,27 W 55* 3 /.2 2346+33.2(54.5) 1.62±0.059 3800±88 

C 2 H 6 

23 v 7 [50] 3.25±0.376(0.556) [1.62] 5.27+0.92 

CH 3 OH f 

CH 3 OH f 

23 v 2 [-50] 30.1+3.07(4.37) r , 48.817.1 

[1.62] 

22 v 3 48 +10 /- 7 28.612.82(3.41) 46.315.5 


a Rotational temperatures for H2O using lines in orders 26 and 27 (Fig. 2), and for 
CH3OH from the V3 band (Fig. 4). Quantities in square brackets denote adopted 
values. All uncertainties in the Table are la . 


b Nucleus-centered production rate based on signal contained in a 9-row extract 
centered on each emission profile peak, as determined by a best-fit Gauss function. 
For H2O (Fig. IB), the emission is centered on the continuum peak (taken to be 
centered on the nucleus), while for C2H6 and CH3OH extracts centered -2 rows (~ 400 
km) from the nucleus produced the maximum emission intensities (Figs. 1C, ID). For 
each tabulated entry, the first uncertainty is the stochastic noise and the second (in 
parentheses) is the standard error among the lines (i.e., spectral intervals) sampled, 

c Growth Factor (GF, Eq. 1) used to correct Q» c for slit losses, which are due primarily to 
atmospheric seeing. For our analysis we assume the value for H2O applies to the 
other molecules. 

d Total or global Q (Eq. 1). For each entry, the uncertainty incorporates the (dominant) 
standard error in Qnc and the error in GF measured for H2O (aGF contributes a 
relatively small amount). Using the GF measured for C2H6 (1.90±0.285, based on the 
profile in Fig. 1C) increases its global Q (and hence the ethane abundance, see Table 
3) by ~17 percent, to (6.17±1.40) x 10 25 S' 1 ; however this is within la of that based 
on the GF from H2O and does not change our conclusion of highly depleted C2H6 in 
comet 21P/GZ. Combining CH3OH from orders 22 and 23 results in Qtot(CHsOH) = 
(47.5+4.4) x 10 25 s' 1 . For the combined CH3OH profile (Fig. ID) we measured 
GF=1. 65+0. 241; using this GF increases the combined Qtot(CHsOH) by only a 
negligible amount (1.9 percent). For CH3OH in order 23, aQ nc and aQ Co t also include 
uncertainties in empirical g-factor summed over each interval, based on stochastic 
noise levels in the 8P/Tuttle residuals (Fig. 3A; see also Table Al). 

e The 31 H2O lines included in our analysis arise from eight distinct non-resonant 
("hot") bands: 101-110, 111-110, 200-001, 210-110, 200-100, 003-002, 101-001, 
and 201-200. 

f For CH3OH in order 23, Q nc and Q to t correspond to the weighted mean of all intervals 
in Fig. 3C excluding number 5, due primarily to (unknown) uncertainties in g-factor 
for the co-measured OH* P12.5 2- line, both for 8P/Tuttle and 21P/GZ as discussed in 
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§4,2.2. Including interval 5 lowers the global Q for CH 3 OH in order 23 by ~18 
percent, to [40.1+6.6] x 10 25 molecules s' 1 , and lowers the mean Q(CH 30 H] (and 
hence the abundance of CH 3 OH from combining orders 22 and 23] by approximately 
10 percent (see Table 3, note f). 
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Table 3. Volatile Abundances in 2 1P/G-Z and Other Jupiter Family Comets 


Comet 

Year 3 

_C2ik a 

h 2 o 

CH^OH 3 
H2O 

Log C2 b 
OH 

C 2 /OH9 

(%) 

C 2 H 6 q 

c 2 

21P/GZ 

1998 

0.22±0.13 c 






1998 

< 0.05-0.08 d 

0.9-1.4 d 

-3.41W 

.056 



1998 


1.6±0.4 e 

-3.14i 2 ) 




2005 

0.139±0.024 f 

1.25±0.12 f 



2.7 

9P/Tempel 1 

2005 

0.23+0.048 

0.55±0.09 h 

1.4+0.28 

1.1+0.28 

-2.901 1 ) 

-3.1712) 

.097 

2.6 

73P/SW3-B 

2006 

<0.3' 





73P/SW3-B 


0.14+0.0H 

0.25±0.04i 

-3.2912) 

.051 

2.7 

73P/SW3-C 


0.15+0.04' 

<0.38‘ 




73P/SW3-C 


0.11±0.01i 

0,21+0.04) 




17P/Holmes 

2007 

1.78±0.26 k 

3.2±0.6 k 

-2.42( 3 ) 

0.38 

5.2 

6P/d'Arrest 

2008 

0.26+0.06 1 

1.99+0. 42 ! 

-2.371D 

-2.9012) 

0.28 

1.0 

103P/Hartley 2 

2010 

0.75±0.03 m 

0.65+0.09" 

1.87±0.13 m 

1.72+0.11" 

-2.361 1 ) 

-2.5912) 


2.3 

lOP/Tempel 2 

2010 

0.39+0.04P 

1.58+0.23P 

-2.701 1 ) 

-2.7112) 

0.20 

2.1 


a Year observed (excluding optical observations; see note b) and abundances of 
ethane and methanol with respect to H2O = 100 (all uncertainties in the Table 
represent la, and upper limits are 3a). Previously published values for 
Q(CH30H) from IR observations using the V3 Q-branch intensity are inter- 
normalized based on g(Q-br) = 1.5 x 10+ ph s' 1 molecule' 1 (Villanueva et al. 
2012 ). 

b Abundances of C2 are taken from the following papers: (1) A’Hearn et al. 1995, (2) 
Fink 2009, and (3) Schleicher 2009. With the exception of 17P, the optical 
observations refer to apparitions different from those listed under "year"; e.g., 
years of observation for Fink et al.: 21P . 1985; 9P, 1994; 73P, 1990 and 1995; 6P, 
1995; 103P. 1997 and 1998; 10P, 1988. A’Hearn et al. 1995 classify comets 
having -1.22 < log(C 2 /CN) < -0.21, and -4.13 < log(C 2 /OH) < -2.98 as "carbon- 
chain depleted", and those having -0.09 < log(C2/CN) < 0.29 and -2.90 < 
log(C2/OH) < -2.10 as "typical". The values for 17P refer to the mean of 2007 Nov. 
01.16, Nov. 20.12, and 2008 Jan. 01.23, during the decline from its spectacular 
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outburst in late October 2007. Its rapidly evolving coma complicates a 
comparison of abundances for primary (parent] and product species, and in this 
sense 17P represents a special case. 

c Mumma et al. 2000. 

d Weaver et al. 1999. Their value for Q(CH 30 H] was reported as a 3a upper limit. It 
was based on the entire V 3 band using g(v 3 ) = 1.5 x 10 ' 4 s' 1 , therefore no re-scaling 
of production rate from the published value is required. 

e Biver et al. 2002. 

f This work. The value listed for C 2 H 6 is based on six V 7 Q-branches (see Figs. 3B 
and 3D and related discussion in §4.2.3]. Including the seventh Q-branch ( r Q3] 
modifies its abundance to (0.136 ± 0.021] x 10' 2 . The value listed for CH 3 OH is 
based on signal in 36 spectral intervals in order 22 (Fig. 4] and 17 spectral 
intervals (i.e., with number 5 omitted] in order 23 (Figs. 3C and 3E). Including 
interval 5 in order 23 results in a mean CH 3 OH abundance of (1.12 ± 0.11] x 10' 2 . 

g Mumma et al. 2005. The first row corresponds to pre-impact, the second row to 
the D1 ejecta. 

h DiSanti et al. 2007b. 

' Villanueva et al. 2006. 

i Dello Russo et al. 2007 and Kobayashi et al. 2007 reported consistent values for 
fragment B, from independent measurements. 

k Dello Russo et al. 2008. 

I Dello Russo et al. 2009. 

m Mumma et al. 2011. 

II Dello Russo et al. 2011. 

p Paganini et al. 2012. 

* For comets with multiple entries for log(C 2 /OH], the quantity C 2 /OH refers to 
their mean value. In calculating C 2 H 6 /C 2 , the relation Q(H20] = 1.1 x Q(OH] is 
assumed. The entry for 21 P is based on the value for C 2 H 6 /H 2 O reported here 
(0.139 percent]. The entry for 9P is based on the pre-impact C 2 H 6 /H 2 O (note g). 
The entry for 73P uses the weighted mean value from the three measurements of 
C 2 H 6 /H 2 O having la errors (~ 0.14 percent]. The entry for 103P uses the 
weighted mean value from the two measurements listed for C 2 H 6 /H 2 O. 
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Figure captions 

Fig. 1. A. Schematics showing the observing geometry for 21P/GZ on UT 2005 June 03. 
Upper : PA of the extended Sun-comet vector (260°) and slit orientation (157°, or 23° W 
of N) projected onto the sky plane, with the comet located at the origin. Lower : Solar 
phase angle ((3 = 44°) for 2 1 P/G-Z. B-D. Comparison of spatial profiles in Comet 2 1 P/G- 
Z. Projected distance is reckoned with respect to the peak of the observed continuum 
profile. For comparison, a representative point-spread-function (PSF) is also shown. 

Fig. 2. Spectral extracts of H 2 0 in Comet 21 P/G-Z, with best fit fluorescence model, 
indicating T rol = 55 +3 /. 2 K and OPR = 3.00+0.15. The prominent feature near 3485 cm ' 1 
(marked with asterisk) corresponds to two blended OH lines (v=l-0: P2.5 f and P2.5 1 + , 
at 3484.59 and 3484.74 cm' 1 , respectively). 

Fig. 3. A. Top traces: Continuum-subtracted residuals from order 23 in 8 P/Tuttle, along 
with our optimized empirical fluorescence model for CH 3 OH. Subtracting modeled 
emissions for C 2 Pl 6 , OH, and CH 4 leaves primarily v 2 emissions (hence labeled “CH 3 OH 
residuals”). Regions marked V contain excess intensity, perhaps due to unaccounted-for 
CH 3 OH V9 emissions, and so are omitted from our v 2 analysis (see Appendix). A stick 
spectrum for all lines included in the model is also shown, both here and for 21P/GZ, 
each being multiplied by their appropriate monochromatic transmittances for direct 
comparison with observed emissions. B. Top trace: Extracted spectrum in 21P/GZ, with 
modeled transmittance function convolved to \/A\ = 24,000. Lower traces: Continuum- 
subtracted spectrum and convolved quantum band model for C 2 H 6 V 7 (with seven Q- 
branches indicated) and empirical model for OH prompt emission. C. Residual emission 
spectrum after subtraction of modeled C 2 H 6 and OH, with our new empirical model 
super-imposed. D. Production rates derived from individual C 2 H 6 Q-branches at 
(adopted) T ro i of 50 K, and the mean Q to t based on points 2 through 6 (dashed line). E. 
Similar graphic for the 1 8 spectral intervals used in our analysis of CH 3 OH in order 23 
from panel C. The mean Q lol excludes point number 5, as discussed in the text. 

Fig. 4. A. Spectral extracts of order 22, showing emissions from the V 3 band of CH 3 OH 
in Comet 21 P/G-Z, with the Q-branch indicated. B. Measured (total) production rates 
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derived from each of the 36 spectral intervals delineated in panel A, with that from the Q- 
branch indicated (Qbr)- C. Best-fit slope over a range of rotational temperatures, 
indicating optimal T ro t (slope = 0) and ±la uncertainties (these are based on 
corresponding uncertainties in slope; cf. DiSanti et al. 2006 for detailed formalism). 
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Fig. 1, A. Schematics showing the 
observing geometry for 21 P/GZ on 
UT 2005 June 03. Upper : PA of the 
extended Sun-comet vector (260°) 
and slit orientation (157°, or 23° W 
of N) projected onto the sky plane, 
with the comet located at the origin. 
Lower : Solar phase angle ((3 = 44°) 
for 21P/G-Z. B-D. Comparison of 
spatial profiles in Comet 21 P/G-Z. 
Projected distance is reckoned with 
respect to the peak of the observed 
continuum profile. For comparison, 
a representative point-spread- 
function (PSF) is also shown. 
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Fig. 2. Spectral extracts of H 2 0 in Comet 21P/G-Z, with best fit fluorescence model, indicating 
T r0 i = 5 5 4 3 /.2 K and OPR = 3.00±0.15. The prominent feature near 3485 cm' 1 (marked with 
asterisk) corresponds to two blended OH lines (v=l-0: P2.5 1* and P2.5 1 + , at 3484.59 and 3484,74 
cm' 1 , respectively). 
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Fig. 3. A. Top traces: Continuum-subtracted residuals from order 23 in 8P/Tuttle, along with our 
optimized empirical fluorescence model for CH 3 OH. Subtracting modeled emissions for C 2 H 6) 
OH, and CH 4 leaves primarily v 2 emissions (hence labeled “CH 3 OH residuals”). Regions 
marked ‘x’ contain excess intensity, perhaps due to unaccounted-for CM 3 OH v 9 emissions, and 
so are omitted from our v 2 analysis (see Appendix). A stick spectrum for all lines included in the 
model is also shown, both here and for 21P/GZ, each being multiplied by their appropriate 
monochromatic transmittances for direct comparison with observed emissions. B. Top trace: 
Extracted spectrum in 21P/GZ, with modeled transmittance function convolved to X/AX - 
24,000. Lower traces: Continuum-subtracted spectrum and convolved quantum band model for 
C 2 H 6 v 7 (with seven Q-branches indicated) and empirical model for OH prompt emission. C. 
Residual emission spectrum after subtraction of modeled C 2 H 6 and OH, with our new empirical 
model super-imposed. D. Production rates derived from individual C 2 H 6 Q-branches at (adopted) 
T [0l of 50 K, and the mean Q, ol based on points 2 through 6 (dashed line). E. Similar graphic for 
the 18 spectral intervals used in our analysis of CH 3 OH in order 23 from pane! C. The mean Q 10l 
excludes point number 5, as discussed in the text. 
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Fig. 4. A. Spectral extracts of order 22, showing emissions from the 
V 3 band of CH3OH in Comet 21P/G-Z, with the Q-branch indicated. 
B. Measured (total) production rates derived from each of the 36 
spectral intervals delineated in panel A, with that from the Q-branch 
indicated (Qb,). C. Best-fit slope over a range of rotational 
temperatures, indicating optimal T rot (slope = 0) and ±1 a uncertainties 
(these are based on corresponding uncertainties in slope; cf. DiSanti 
et al. 2006 for detailed formalism). 
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Appendix. Application of the CH3OH Empirical v 2 Band Model to Comet GZ 

For our empirical model of the CH3OH v 2 band, we adopt notation according to the G6 
molecular symmetry group, as this lakes into account the large amplitude torsional 
tunneling between the three equivalent CH 3 orientations relative to OH (Xu et al. 1997; 
Xu et al. 2008). Under the G6 group representation, the rotational structure is described 
using three quantum numbers (J, K, /), where J (> 0) is the usual total angular 
momentum, K is its projection along the symmetry a-axis of the molecule (such that 0 
< K < J), and / is related to the parity of the states involved in the transition. The G6 
symmetry group uses the representations A], A 2 , E| and E 2 , where subscript “1” is 
associated with positive parity of a ro-vibrational slate and “2” with negative parity, 
thereby leading to the representations A' h , A', E + , and E'. Following Mekhticv et al. 
(1999), we differentiate between E" 1 to E' states by denoting the latter using negative K 
values. Using the standard convention, a single-prime denotes the upper vibrational state 
(v = v 2 ) and a double-prime denotes the lower (ground) vibrational state (v - 0). Then, 
for example (see Table A. 1), the designation A 1 (1->0)P4 corresponds to J' = 3, K r = 1, /' 
- 1, and J" = 4, K" = 0, /" = 1, while E (-l->-2)R5 corresponds to J’ = 6, K' = -1, /’ ~ 2, 
and J" - 5, K" = -2, l" - 2. 

To establish spectral intervals for determining the CH3OH production rate using our 
empirical model, we first identified regions in the order 23 spectrum of Comet 8P/Tuttle 
having residual intensity after accounting for emissions due to C 2 H<s, OH, and CII4 (see 
Fig. 3A, in particular the bottom trace labeled “CH3OH residuals”). We next 
incorporated line assignments and frequencies from the compilation in Xu et al. (1997) 
based on supersonic (jet) spectra of CH 3 OH measured at low temperature (~ 17 K). We 
assigned fluorescence g-factors for lines contained in each spectral interval such that their 
sum reproduced the nucleus-centered production rate (Q nCl /«/ in Eq. 2) retrieved for 
8P/Tuttle (=2.90 x 10 26 molecules s' 1 ) based on application of our recently published v 3 
quantum band model to the residuals shown in Fig. 4A for order 22 (this “revised” Q,ot, 
which includes GF = 1.5 as discussed in Bonev et al. 2008, is presented in Table 4 of 
Villanueva et al. 2012; also see §3.2 in the present paper). 



Rigorous calculation of line g- factors requires a full quantum mechanical fluorescence 
model for the V 2 band, in which solar pumps from all (dipole-allowed) contributing levels 
in the ground vibrational state (v = 0) into each upper level (in v = V 2 ) are summed (with 
potential cascades from higher ro-vibrational levels included). Each particular line g- 
factor is then given by this summed pump rate multiplied by the rotational branching 
ratio for the line in question. However, such a treatment is beyond the scope of the 
present paper. 

Instead, in our empirical treatment, we preformed checks on lines having similar 
rotational designation (J, K) and symmetry (A, E). This approach provided more 
consistent g-factors, and also provided insights as to which spectral regions not to include 
in our treatment. To accomplish this, we examined the intensities of spectral lines, which 
for a given transition can be expressed as a proportionality with respect to line parameters 
and rotational temperature: 

S Hne o : ^-L hl [1 - exp(-hcv/kT)][(2J +l)exp(-hcE/kT)] 

Here, and in the formulae below, v is the line frequency (cm' 1 ), vo is the band frequency 
(2999 cm' 1 ), J, K and E (cm' 1 ) are respectively lower state rotational quantum numbers 
and energy (denoted “Ei ow ” in the table below, and taken from Mekhtiev et al. 1999), T is 
rotational temperature, and Lhl is the Honl-London factor. For a perpendicular band 
such as V 2 , the intensity of which is dominated by lines having ‘A 5 symmetry in b-type 
sub-bands (for which AK - ±1 ; e.g., see Xu et al. 1997), Lhl is given by: 

r (J-l + K)(J + K) AT w _ , ,. N 

L hl = — for AJ = -1 (P-branch lines), 

J *4" J ) 

L hi = - — ■■■ “ + ^ for AJ = 0 (Q-branch lines), (A2) 

J ( J *"H 1 ) 

T (J +2 ± K)(J + 1 ± K) , 1T , _ , , . 

L hi = for AJ = +1 (R-branch lines), 

,,L (J + 1)(2J + 1) 

where the upper sign in each equation refers to AK (= K upper -> K /oi( , c ,.) = +1, and the 
lower sign refers to AK = -1 (Herzberg 1945, p. 426). 



An example is provided by P-branch lines in the A(AK = 1 -> 2) sub-band. In the Comet 
8 P/Tuttle “CH 3 OH residual” spectrum (Fig. 3A), the intensity contained in the left-most 
(highest frequency) peak in interval 16 (near 2976.0 cm' 1 ) is dominated by A' /+ ' > ' /+ P4 
components. The corresponding P3 components at ~ 2977.6 cm ' 1 coincide approximately 
with a feature having strong residual intensity (additionally, a second feature near 2977.4 
cm ' 1 has lower intensity but no viable corresponding V 2 line). Using the above equations 
(and assuming T = 50 K) suggests an intensity ratio P3/P4 = 1.07, and this is reflected in 
the model shown in Fig. 3. [We note that the modeled “P3” feature appears weaker than 
the feature in interval 1 6, however this results from its relatively poorer transmittance (~ 
50 %) compared with that of the “P4” feature (~ 95 %).] However, to account for the 
intensity in the CH3OH residuals near 2977.6 cm ’ 1 requires an unrealistically large ratio 
of ~ 1.8, and we therefore exclude this region from our analysis. Similar reasoning 
pertains to other excluded spectral regions (marked ‘x’ in Fig. 3A), and this approach will 
be used as our empirical model is improved through extension to other comets in our 
database encompassing a range of rotational temperatures. Table A1 provides a summary 
of our empirical model as applied to Comet 21P/Giacobini-Zinner. 



Table Al. CH3OH V 2 empirical model results for 21P/Giacobini-Zinner on UT 2005 June 03 


Or Vma X -V mi n* Rot ID b 


(cm' 1 ) 


Restv b Eiow 


wtd mean 


g-fac d Qe 


CIO- 7 s 1 } (1025 s- 1 ) 


A 1 (0->0) P3 3002.154 


A* (l->0) P4 

3002.121 

E(l->2) R6 

3002.054 

E(l->0] R2 

3001.975 

A**>*(3->3)Q 

3001.940 ; 


3001.905 

E[-3] R6 

3001.799 

A-(1->1) R1 

3001.775 

A* (0->l) R2 

3001.777 

| A' (1->1j R1 

3001.735 



E(-1->0)Q 

3001.080 

E(0-0]R3 

1 3001.042 


9.68 


16.13 


49.35 


13.96 


-50.00 


42.21 


11.73 


14.90 


77.47 


11.70 


3001.954 25.65 0.831±0.284 15.72*1.91 11.10*3.79(4.04) 


3001.079 37.65 0.790*0.210 32.69*1.51 50.73*13.49(13.82) 




A*/-(2->2)R5 


E(-l-> -1)P4' 


E(-l-> -1)P4" 


E(0-> -1)Q 


3000.002 


2999.912 


2999.894 


2999.836 


2999.83 


2999.854 24.31 0.625*0.202 23.36*1.36 56.19*18.21(18.61) 




E(1->2)R3 

2997.250 

A* (0-l)P7 

2997.161 

E(1->0)Q 

2997.14 

E(0->1)R4 

2997.130 

E(-2->-2}P3 

2997.116 


2997.139 33.02 1.355*0.191 37.21*1.30 76.59*10.82(11.49) 



























































































A 4 • > ’{1->0}Q2 

2996.870 

A*->-(l->0)Ql 

2996.847 

A 4 ->-(l->0)Q3 

2996.811 

A 4 ->-(l->0]Q4 

2996.763 

E(-l-> -1)P6" 

2996.715 

A 4 ->-(l->0)Q5 

2996.702 

A 4 /--> 4 /-(2->2)R3 

2996.674 

E(0--> -1)P2‘ 

2996.639 


E(-l-> -1]P6' 


E(0-> -1)P2" 


2996.626 


2996.611 


2996.542 


A’->-(l->0)Q7 2996.537 


E(1->1)R2 


2996.470 


A*->-(0->l)Q8 2996.431 


2996.675 32.01 0.809±0.289 85.88^=1.98 19.82*7.08(7.13) 


A - (1— >2}R5 

2992.389 

E(0-1]R1 

2992.324 

E(1-0)P3 

2992.294 

E(-l-> -2)Q 

2992.23 

A 4 (1->1}P4 

2992.153 

A 4 (0->l]P3 

2992.131 


2992.224 45.85 0.689*0.409 45.83*2.67 31.68*18.81(18.94) 


2990.89 - 2990.53 E(1->2)Q 2990.76 


A-(1->2)R4 2990.693 


E (1->0)P4 2990.684 


A* (1->1)P5 2990.580 


A 4 (0->l)P4 2990.575 


2990.715 42.30 0.511*0.246 36.60*1.67 29.43*14.19(14.29) 


A 4 (0->l)P7 

2986.017 

E(-2-> -3)R5 

2985.956 

E(1->2)P3 

2985.956 

E(0->1)P2‘ 

2985.887 











































































































E( 1 -> 0) P 7 

2985.859 

E(0->1)P2" 

2985.859 

E (’!•> -2)P4 

2985.708 


2985.907 


0.200±0.279 23.36^1.90 18.10±25. 19(25.24) 


2984.62-2984.33 A*(0->1)P8 


E(1->2)P4 


2984.524 


2984.382 


2984.342 


2984.403 50.99 0.157±0.317 13.61±1.76 24.32±49.23(49.33) 




A~>*(1->2)Q10 

2983.131 

A-->’(1->2)Q9 : 

2983.010 

(0 — > 1 )P9 

2983.044 

A" >+ (1->2)Q8 

2982.898 


2983.019 92.44 0.755±0.206 17.91±1.43 89.14±24.3 1(25.54) 




A + -> • ( l-> 2} Q7 

2982.797 

E(-2-> -3)R3 

2982.797 

E(1->2)P5 

2982.726 

A’ ■> ■ (1->2)Q6 

2982.705 

E(0->1)P4 

2982.635 

E{1->0)P9 

2982.635 

A->-(l->2)Q5 : 

2982.624 

A* ■>' (1->2)Q4 : 

2982.555 

E(-l-> -2)P6‘ 

2982.573 

E(-l-> -2)P6" 

2982.487 

A‘->-(l->2)Q3 | 

2982.501 


2982.647 49.90 0.444±0.320 50.57±1.89 18.55±13.39(13.43) 




E(1-1)P6 

2981.941 

A + ->-(3->2)Q 

2981.86 

A--> + (3->2)Q 

2981.83 


2981.848 78.03 0.683±0.206 39.17±1.40 36.89±1 1.14(11.30) 


2981.21 -2980.76 E(1->2)P6 


E(0->1)P5 


2981.106 


2981.012 


2981.04 












































































































E(-l-> -2JP7 

2980.924 

A” >4 (3->4)Q 

2980.825 

A->-(3->4)Q 

2980.793 



E(0-> -1)P12 

2979.549 

E(1->2)P7 

2979.479 

E(0->1)P6 

2979.380 

E(-l-> -2)P8 

2979.332 

A 4 (1->2)P2 

2979.181 

A-(1->2)P2 

2979.204 


2980.918 60.21 1.512±0 .666 44.96*4.91 71.10*31.33(32.38) 


2979.352 54.16 1.638*0.641 24.33*3.06 142.40*58.34(61.25) 


2977.16-2976.92 E(1->1)P9 


A-(2->3)P3 


2977.044 87.24 


2976.965 44.29 


2976.966 44.29 


2976.973 46.42 0.637*0.224 11.36*1.58 118.70*41.78(45.13) 




E(-l-> -2JP10 

2976,192 

E(0->1)P8 

2976.043 

A-(1->2)P4 

2976.047 

A 4 (1->2)P4 

2975.909 

A 4 (2->2)P9 

2975.750 

A- (2->2)P9 

2975.735 

E(2->3)P7 

2975.656 

A-(2->3)P4 

2975.356 

A 4 (2->3)P4 

2975.357 


2975.732 57.00 1.800*0.422 36.29*2.81 104.98*24.61(26.20) 


2974.35 -2974.14 E (0->l)P9 


A 4 (1->2)P4 


A-(3->4)P4 


2974.282 


2974.260 


2974.240 


2974.239 




A- (2->3)P5 

2973.748 

A 4 (2->3)P5 

2973.751 


2974.251 67.48 -0.079*0.222 16.63*1.95 -10.01*28.30(28.33) 


58.81 




















































































































2973.750 

58.81 

0.544±0.275 

8.70*1.70 

132.58*67.07(72.06) 



— 









a Range of rest frequencies (cm' 1 ) encompassed by each spectral interval in Fig. 3. 

b Rotational designation and tabulated rest frequencies are taken from Xu et al. 1997. The asterisk in 
interval 13 denotes an H-type Q-branch of the vg band; its frequency and energy are only 
approximate. The heading “wtd mean” for each spectral interval denotes the mean rest frequency 
weighted by contributing line-by-line empirical g-factors (including atmospheric transmittance at the 
Doppler shifted line center, using A-dot = -5.42 km s' 1 ), based on the “CH3OH residuals” in 
8P/Tuttle (bottom trace in Fig. 3A). Equations SOM 1 and 2 were used as a guide for modeling line 
intensities for lines having similar rotational designation and symmetry, as discussed above. Lower 
state energies (E| 0W ) were taken from Mekhtiev et al. (1999). 

c Nucleus-centered line flux and la stochastic error for each interval, corrected for weighted mean 
atmospheric transmittance. 

(i Summed fluorescence-efficiency g-factor (at R h = 1 AU) for each interval. Uncertainties in g 
represent la and are based on the stochastic noise summed over for each spectral interval in the 
8P/Tuttle order 23 spectrum. 

e Global (total) production rate based on the signal contained in each interval. The first uncertainty 
incorporates the la error in line flux (plus aGF, which adds only a minor contribution). The second 
uncertainty (in parentheses) includes these errors plus the la uncertainty in empirical g-factor for 
each interval. 
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